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ABSTRACT 

We present self-consistent calculations of non-equilibrium (time-dependent) cooling 
rates for a dust-free collisionally controlled gas in wide temperature (10 K ^ T ^ 
10^ K) and metallicity (lO""^ Zq Z 2 Zq) ranges. We confirm that molecular 
hydrogen dominates cooling at 10^ < T < 10" K and Z < lO'^ Zq. We find that the 
contribution from H2 into cooling rate around T ~ (4 — 5) x 10"^ K stimulates thermal 
instability in the metallicity range Z ^ 10^^ Zq. Isobaric cooling rates are gener- 
ally lower than isochoric ones, because the associated increase of gas density leads 
to both more efficient hydrogen recombination and equilibration of the fine-structure 
level populations. Isochoric cooling keeps the ionization fraction remains quite high at 
T < 10"* K: up to ~ 0.01 at T ~ 10^ K and Z < 0.1 Zq, and even higher at higher 
metallicity, contrary to isobaric cooling where it at least an order of magnitude lower. 
Despite this increase in ionization fraction the gas-phase formation rate of molecu- 
lar hydrogen (via H~) lowers with metallicity, because higher metallicity shorttens 
the evolution time. We implement our self-consistent cooling rates into the multi- 
dimensional parallel code ZEUS-MP in order to simulate evolution of a supernova 
remnant, and compare it with an analogous model with tabulated cooling rates pub- 
lished in previous works. We find significant differences between the two descriptions, 
which may appear, e.g., in mixing of the ejected metals in the circumstellar medium. 

Key words: atomic processes - molecular processes - plasmas - galaxies: general - 
intergalactic medium - ISM 



1 INTRODUCTION 

A choice between equilibrium and non-equilibrium 
(time-dependent) cooling rates in collisionally controlled 
gas strongly depends on the ratio between dynamical and 
chemical/cooling times. Calculations of the cooling rates 
of astrophysical plasma in the coUisional ionization equi- 
librium (CIE) for T ^ lO'^ K were performed by many au- 



thors jHousd Il964l: ICox fc Tucked [l969: Raymond et al. I 
19761: IShuU fc van SteenberelllQSd: [Caetz fc Salpetd | l983| : 



Bohringer fc Hensled 19891: ISutherland fc Dopital Il993l: 
Landi fc Landinill 19991: iBeniamin et al. ll200ll :l Brvans et al. I 



20061 : ISchure et al. I 120091 ^ However, calculations of the 



time-dependent ionization of metals and associated ra- 



the CIE states and rates (IKafatod 19731: Shaoiro fc Moore 


1976 


: lEdsar fc Chevalieij 1 19861: ISchmutzler fc Tscharnutei 


1993 


; 'Sutherland fc Dopita' 


1993': 'Gnat fc Sternberg 20071: 


de Avillcz fc Brcitschwcrdt..20ia:iVasilievi,2011). 



A more complex task is a calculation of cooling 
rate below 10* K, because we need to take into ac- 
count iMJij;_2rocesses_betwe^^ 

(e.g. iDalgarno fc McCra^ll972l : iHoUenbach fc McKedll979l . 
Il989h . Apart from these difficulties a level of the calculation 
complexity of the cooling rate depends on one's purposes 
and typical timescales of physical problem considered. Under 
some conditions chemical network a t T < 10* K can be con- 
strained signi f icantl y. For example, lAbel et al. I (|l997l ) and 
ICalli fc Fallal (|l998l ) pr oposed the reduced chem ical network 
for the primordial gas, IClover fc JappsenI (|2007l ') considered 
the chemistry of the metal- enriched interstellar medium. 
Sometimes, only cooling rates (without solving chemical ki- 
netics equations) are used to study the interstellar medium, 
for in stance, the rates calculated by IDalgarno fc McCravl 
(|l972l ) are applied widely. 

Usually the cooling rates are provided separately 
in two tem perature ranges, the fi rst is T = 10* — 
10* K fe.g.. iRavmond et al. I Il976l: ISutherland fc Dopital 



I1993I : iBrvans et al. I I2OO6I : ICnat fc Sternberg! 12007]) and 
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the second is T 10 K (e.g., IDalgarno fc McCrav 



[l97i; 
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iHoUenbach fc McKed [l989lV Very rarely the rates are pre- 



Spaans fc Norman|[l997l ') . This is 



Table 1. Solar elemental abundances 



sented in both ranges (e.g. 
explained by difference in physical problems considered by 
authors. But very frequently one needs cooling rates in both 
ranges. So that usually high- and low-temperature cooling 
rates are combined from different tables, where the rates 
were computed using different assumptions. Thus, it is nec- 
essary to conduct a self-consistent calculation of cooling rate 
within both high- and low-temperature ranges. 

In this paper we present self-consistent calculations of 
non-equilibrium cooling rates of a coUisionally cooled metal- 
enriched dust-free gas in wide temperature range, fO ^5 
T 10* K, and we study dependence of the rates on metal- 
licity and density. A possible influence from dust grains and 
external ionizing radiation field will be considered elsewhere. 
The paper is organized as follows. In Section 2 we briefly de- 
scribe the details of the model. In Sections 3 and 4 we present 
our results. In Section 5 we study the supernova shell evo- 
lution using newly calculated cooling rates. In Section 6 we 
summarize our results. 



2 MODEL DESCRIPTION 

The chemical and thermal evolution of a gas parcel can 
be divided into a high-temperature (T > 2 x 10* K) and 
a low-temperature (T ^ 2 x 10** K) ranges. Such division 
can be explained by the transition to neutral gas and the 
formation of molecules in the latter range. The full descrip- 
tion of the method and the references to the ato mic data 
for th e former temperature range can be found in IVasilievI 
l|201lf ). Briefly, in our calculations we consider all ionization 
states of the elements H, He, C, N, O, Ne, Mg, Si and Fe. 
We take into account the following major processes in a col- 
lisional gas: coUisional ionization, radiative and dielectronic 
recombination as well as charge transfer in collisions with 
hydrogen and helium atoms and ions. 

To study the chemical evolution of a low-temperature 
(T < 2 X 10"* K) gas we should take into account molecules, 
especially molecular hydrogen, which is expected to be the 
most significant coolant in a gas with metallicity below 
1O~''Z0. Hence, the above-listed ionization states of the 
elements is replenished by a standard set of species: H~, 
H2, HJ, D, D+ D~, hp, needed to model the H2/HD gas- 
phase kinetics l|Abel et al. Ill997l : iGaUi fc Fallal Il998|). Th e 
corresponding rates are taken from iGaUi fc Fallal (|l998h : 
IStancil et aU (|l998l ): iMac Low fc Shulll l| 19861 ). The hst of 
the chemical reactions for molecular hydrogen and deu- 
terium chemistry is presented in Appendix 1. 

The system of time-dependent ionization state equa- 
tions should be complemented by the energy equation, 
which can be written in the form 



du 
H 



p_dp 
/)2 dt 



L 



(1) 



where u = p/(7 — l)p is the thermal energy, p is the gas 
pressure, p is the gas density, 7 is the adiabatic index, and 
L is the generalized energy loss function. The equation {T]) 
can be expressed as 



dp 



1 



7 — Irft J — 1 p dt 



^dp 



(2) 



Element 



(X/H) 







Helium 
Carbon 
Nitrogen 
Oxygen 
Neon 

Magnesium 

Silicon 

Iron 



0.081 
2.45 X 10-4 
6.03 X 10--"^ 
4.57 X 10-4 
1.95 X 10-4 
3.39 X 10-5 
3.24 X 10-5 
2.82 X 10-5 



Assumin g 7 = 5/3 we fi nd that the gas temperature evolves 
as (e.g., lKafatos|[l973 ) 



dT 
'dt 



rieriff A T dn 
Anks n dt 



(3) 



where n is the total particle density of the gas, ne and 
uh are the electron and total hydrogen number densities, 
A{xi,T, Z) is cooling rate, fcs is the Boltzman constant, A 
is a constant equal to 3/2 for isochoric and 5/2 for isobaric 
cooling. The second term in equation ((3]l refiects the relative 
change in the number of particles in the unit of volume, such 
change is very important in gas below T ^ 10* K, where gas 
recombines rapidly. 

The total (except molecular deposit) cooling rate is 
calcul ated using the photo ionization code CLOUDY (ver. 
08.00, jPerland et anil998l]. Th e radiative losses by H2 is 
taken from ( Galli fc Fallal 19981) and the HP co oling rate 
is taken from (|Floweil I2OO0I : iLipovka et al.1l2005h . For the 
solar metallicity we adopt the abundances reported by 
lAsplund et al. I (|2005l ). except Ne for which t he enhanced 
abundance is adopted jPrake fc Testal l2005h . the solar 
abundances are listed in Table [T] In our calculations we 
assume the helium mass fraction Yho ~ 0.24, which cor- 
resp onds to [He/H]=0.081 . Such value is close to observed 
one (jlzotov fc Thuanll 19981 ). In general, it would be difficult 
at this time to ex clude any v alue of Yho inside the range 
0.232 - 0.258 (e.g. IOlivell2010l ). We solve a set of 103 cou- 
pled equations (102 for ions, atoms and molecules, and one 
for temperature) using a Variable-coefficient O rdinary Pif- 
ferential Equation solveiQ (|Brown et al. |[l989l ). 

In our calculations we don't include the chemical ki- 
netics of carbon monoxide and lithium hydrid, dust and 
photoionization effects. It is well known that CO and 
H2O molecules significantly affect on the thermal state 
of the interstellar gas fe.g.. piollenbach fc McKed Il979l : 
lOmukai et ani2005l : lGiover fc Jappsenll2007l ). The CO and 
H2O cooling rates dominate over the C a nd O fine-structure 
coolin g in gas with T > 500 K (e.g., [ci over fc JappsenI 
l2007h . but whether it is possible to produce enough CO and 
H2O molecules in the gas during its cooling timescale. In the 
next section we show that the CO/H2O formation timescale 
is greater than the cooling time, so that such molecules can- 
not form in significant amounts within timescale considered 
here. Hence, we can neglect CO/H2O kinetics in this model. 
The LiH molecule is a significant coolant only at T ^: 50 K. 



^ The numerical code 
jhtt p : / / www . netlib . org/ode/ vodofj 



available 



from 
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Figure 1. The cooling rates for solar metallicity. The isochoric 
rate for gas with n = 1 cm~^ calculated in this work is depicted 
by solid line. The other lines correspond to the rates obtained in 
the previous calculations (see text for details). 



More important influence on the thermal and chemical evo- 
lution is expected from dust and ionizing radiation. Dust 
particles confine a significant part of metals, produce addi- 
tional cooling or heating, catalyze the H2 formation as well. 
The ionizing radiation ionizes gas and dust. We are going 
to consider these effects elsewhere. But we note that a sim- 
ple estimate gives that the H2 formation on dust grains is 
inefficient within timescale considered here (see the next sec- 
tion). Several effects from th e ionizing radi ation on gas at 



T > 10* K were described in (|Vasiliev 



^ radi a 
Mil). 



3 ISOCHORIC COOLING RATES 

For isochoric case we start our calculations from T = 
10* K and follow it until the temperature becomes low 
as T = 10 K. In general, the evolution of gas de- 
pends on initi al ionic composition, temperature and den - 
sity (see e.g. [Sutherland fc Dopital 1 19931 : IVasihevI |2012| ). 
but here we will consider cooling of gas from very high 
temperature, so that there i s no depende nce on initial 
ionic composition of gas (e.g. IVasilievI I2OIII . see the next 
section also). Initially the ionization composition corre- 
sponds to the CIE at T = 10** K obtained from CLOUDY. 
Figure [T] presents the isochoric cooling rates for a gas 
with number density n = 1 cm~'^ and solar metallicity. 
The isochoric rate calculated in this work is depicted by 
solid line. To compare our calculations of the rates to 
the previous ones we show the cooling rates obtained by 
[Sutherland fc Dopital (Il993l) (SD93, we have choosen their 
non-equilibrium d at a), ISpaans fc NormanI ( 19971) ( SN97) , 
iGnat fc Sternberd (l2007l) (GS 7). ISchure et al. I (|2009l ) 
(S09). and lDalgarno fc McCravl (|l97a )" for ionization frac- 



tion, fi = n^/uH = 10"", 10"^, 10"^, 0.1, which are 
shown by dash-dot-dotted lines from bottom to top. The 
data obtained from CLOUDY code is depicted by crosses 
(the H2 molecules are ignored in the equilibrium calcu- 
lat ion). One should not e that the cooling rates obtaine d 
by [Sutherland fc Dopital (|l993l ). ISpaans fc NormanI (|l997l ). 
iGnat fc Sternberg! (|2007f ) were got in non-equ librium calcu- 
lations, so that these rates are more or l ess close to our rate 
within temperature range 10"* - 10^ K. ISpaans fc NormanI 




Figure 2. Upper panel. The isochoric cooling rates for gas with 
10"* Zq, 10-3 Zq, 10-2 Zq,0.1 Zq,1 Zq,2 Zq metallicities. 
The gas density is n = 1 cm~^ . Middle panel. The temperature 
dependence on the value of fluence, V = J ^dt (for isochoric pro- 
cess it is simply 77 = tn), for the same gas parameters. Lower 
panel. The isochoric cooling time. 



l|l997h didn't include high ionization states of metals in 
their model, so their cooling rate demonstrates significant 
difference at T > 10*^ K. The deviation at T below lO'" K 
can be explained by different atomic data used in the cal- 
culations, and, as a consequence, difference of the ion- 
ization states of m etals. The cooling curves obtained by 
ISchure et ai~l (|2009l ) and by using the CLOUDY code are 
examples of CIE rates, so that significant difference of these 
rates from non-equilibrium calc ulations is not a surprise. 
The iDalgarno fc McCravl (|l972h cooling curves were ob- 
tained by simply summing up cooling rates of hydrogen 
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log T, K log T, K 



Figure 3. The contributions to the total isochoric coohng rate 
from each chemical element for 10~^ Zq (upper panel) and Zq 
(lower panel). The contribution by molecules (H2 and HD) is 
depicted by grey solid line, the total rate is presented by the 
thickest solid line, the deposit from H and He is shown by thin 
solid line. 



and metals, in which the fractional ionization is used as 
a free para meter. One can note that our cooling rate is 
close to the iDalgarno fc McCravl (|l972l ) rate for Xe = 0.1 
at T ~ 10^ - 10* K and drops to that for = 10"" at 
T ^ 10^ K. This proves a necessity of self-consistent calcu- 
lations of cooling rates. 

Figure [Scupper panel) presents non-equilibrium cooling 
rates for a dust-free coUisionally controlled gas in wide tem- 
perature {IQ K i^T ^ Iff K) and metallicity (10"" Zq < 
Z ^2 Zq) ranges. The middle panel shows the temperature 
evolution (we present the dependence on fluence, = / ndt, 
which is for isochoric process 77 — tn, where t is time elapsed 
from the beginning of the evolution, n is the number density 
of gas). Here we should note that all coUisional processes 
are proportional to the square of density. To see clearly 
the importance of non-equilibrium effects the cooling time, 
tcooi = fc(,r/A(r)n, is depicted at lower panel. One can see 
that the cooling time is comparable with the time elapsed 
since the beginning of the evolution at very high tempera- 
ture, log T ^ 7.5, (early time of the evolution) for any metal- 
licity and at low temperature, log T ^ 2.5 for Z ^ 10~^Zq. 
In other temperature and metallicity ranges the cooling time 
is shorter than the evolution time, thus, the non-equilibrium 
effects are significant. 

Figure [3] shows the contributions to the total cool- 



Figure 4. The electron, Xe, (upper panel) and molecular hy- 
drogen, xjjji (lower panel) fractions for metallicities 10~* Zq, 
10-3 Zq, 10-2 Zq, 0.1 Zq, 1 Zq and 2 Zq in isochorically 
cooled gas. 



ing rate from each chemical element and molecules for 
lO"'^ Zq (upper panel) and IZq (lower panel). Coolants 
at temperature higher tha n 10* K are describe d by e.g. 
ISutherland fc Dopital (119931 ): ISchure et aFI (|2009l ): IVasihevI 
(20111). Also it is well-known that below 10* K carbon, 
oxygen and H2 molecules do minate in cooling depending 
on m etallicity of gas (see, e.g. iHollenbach fc McKed Il979l . 
Il989l ). At very low metallicities, ^ 10 ^ Zq, molecular hy- 
drogen is the main coolant below 10* K. For Z ~ 10"'^ Zq 
the contribution from carbon becomes significant only at 
T ^ 100 K (upper panel). The increase of metallicity leads 
to rise of contribution from metals and shortening of tem- 
perature range, where the H2 dominates in cooling. For 
Z ^ O.l Zq the H2 deposit becomes minor at T < 10* K, 
e.g. for solar metallicity the molecules produce a negligible 
cooling rate, whereas the contribution from metals is domi- 
nant (lower panel). 

One should note that thermal instability can de- 
velop at T < 10* K. For Z < 10"^ Zq the instabili t y cri- 
terion for a gas cooled i sochorically (jShchekinovl Il978l : 
IChevaher fc Iniamuralll982l l 



is fulfilled around a small bump in the cooling rate at 
T ~ (4 — 5) X 10^ K, which is produced due to the contri- 
bution from molecular hydrogen (Figure[3K). The instability 
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criterion is not reached within 10 '^Zq ^ Z ^ Zq. Again the 
conditions favoured to thermal instability can be found for 
Z > at T ~ (0.8 - 7) X 10^ K (Figure Eb), where the 
instability develops due to cooling by metals. 

The higher metallicity of a gas is the faster evolution 
goes, or in other words, the cooling time is decreased with 
the increase of metallicity, e.g. at T ~ (1 - 10) x 10^ K 
the cooling timescale for gas with solar metallicity becomes 
shorter more than ten times compared to that scale for gas 
with Z = 0.1 Zq (lower panel of Figure [21). Protons and 
electrons cannot recombine rapidly and such gas remains 
overionized: the electron fraction is high as 0.03 at lOOK 
for solar metallicity (Figure 2^). Such fast evolution reveals 
in decrease of the H2 fraction (although the fraction is 
high enough and depends on metallicity weakly: it decreases 
from ~ 1.4 X lO"'' for 10"'' Zq to 5 x 10"* for Zq, the H+ 
fraction has the same order), so that for solar metallicity 
the H2 fraction decreases more than 1 00 times compared t o 
its 'universal' value, (2 - 3) x 10"^, (|0h fc Haimanll2002t ). 
which is reached for Z ^ 0.1 Zq (Figure[4j3). 

Here we consider the evolution starting from very 
high temperature, which initially corresponds to the con- 
ditions behind strong shock wave, v > 100 km s"\ Usu- 
ally it is assumed that dust grains existed in a gas be- 
fore such shock wave are destroyed behind it, because 
the destruction time due to thermal sputtering, r ~ 2 x 
10* yr (a/0.01^m)/n (where a is the size of the dust parti- 
cle, see lDraine fc Salpetei|[l979l ). is shorter than the cooling 
time for a gas cooled from T Si 10^ K depicted in Figure [2}5. 
Moreover, the accretion growth r ate is negligible compare d 
to sputtering rate (see Figure 7 in lDraine fc Salpeteilll979l ). 
so that dust grains have no time to reform before the final 
time is reached in our calculations. Hence, in the present cal- 
culations we can do not take into account the H2 formation 
on dust grains. 

But if dust grains nevertheless exist in a gas, then 
it is well-established that they c atalyze the H2 pro- 
ducti on in the interstellar medium (|Hollenbach fc McKej 
1 19791 ). In the warm neutral gas (T > 10^ K) the H2 
formation of grains dominates over that due to the 
gas -phase reactions at Z ^ 0.1 Zq (see the analysis 
m iGloveil I2OO3I '). The time of H2 formation on grain- 
surface is ~ 3 x 10'\T/100)-'^^ inZ/ZQ)-'s 

jHoUenbach fc McKej ioTol fl The H2 formation process can 
be effective, if the formation time is shorter than the time 
elapsed from the beginning of the evolution presented in 
Figure [2}d. More exactly, we are interested in the part of 
the evolution starting from the moment, when the temper- 
ature reaches T ~ 10** K, to the end the calculation, when 
T = 10 K. Just in this temperature range hydrogen recom- 
bines efficiently and H2 molecules can form on grain sur- 
face. This timescale is the same order of magnitude as the 
full time of the calculation (Figure [2}3). One can see that 
the tg^^i^s is about 1-2 orders of magnitude longer than the 



^ A recent anal y sis of H2 formation on grain surfaces by 
ICazaux fc Tielenj ||2004| ) that takes both physisorbed and 
chemisorbed hydrogen into account demonstrates that in the con- 
ditions of interest in this paper the compu ted H2 formation rate 
is ver y similar to the widely used rate of iHoUenbach fc McKej 



time needed to cool from r ~ 10'' K to 10 K. So that, we 
note that neither grain-surface nor gas-phase H2 formation 
channel is effective in a gas with Z > 0.1 Zq, and the H2 
fraction remains negligible (see lower panel of Figure [4]). 

In general, if dust grains survive after strong shock 
waves, they contain a significant part of metals. In this case 
the evolution time of a gas becomes longer and, as a con- 
sequence, grain-surface molecule formation becomes more 
efficient. Because of several additional free parameters, e.g. 
dust destruction efficiency, depletion factor, we will study 
the effects from dust on thermal evolution elsewhere. 

Another molecules, which can significantly affect 
the thermal evolution of gas, are CO and H2 Q (e.g., 
iHoUenbach fc McKee|[l979l : lGlover fc Jappsenll2007l l . These 
molecules mainly form via reaction channels initiate d 
by the hydroxyl OH (e.g., IHoUenbach fc McKeel Il979l ). 
The OH radical forms through reaction between neu- 
tral oxygen and molecular hydrogen. Then, the lack of 
H2 molecules produces smaller OH fraction. Moreover, 
as it is mentioned above the H2 formation on dust 
grains is inefficient in the conditions considered here, and 
its fraction decreases with metallicity. The CO forma- 
tion timescale (for H2O molecule the timescale is close) 
is t(O^CO) ~ (10 "Va:H,) (100cm" Vn)(ZB/Z) Myr 
(|Glover fc Jappsenir2007l 'l. Comparing t{0 -J- CO) with the 
cooling time, which is tcooi ~ 0.01 Myr (Icm^^/n) at 
T ~ 10^ - 10" K for solar metaUicity (Figure [2J;) , it is easy 
to see that the CO and H2O molecules will not form in 
quantities large enough to dominate the cooling. 

The higher gas density is the smaller recombination 
timescales of hydrogen and metals become, so the above- 
described picture is expected to change in isobaric process. 



4 ISOBARIC COOLING RATES 

First of all we should choose initial temperature and den- 
sity of a gas for isobaric calculations. Certainly, to com- 
pare with the isochoric models we should start our calcu- 
lations from T = 10* K. However we extend our calcula- 
tions to T < 10'* K, in this range main radiative losses are 
through fine-structure lines of metals. The level populations 
of these lines go from equilibrium to non-equilibrium within 
n ~ 1 - 100 cm"^ at T <, 10^ K. To catch this transition we 
should set initial density n ^ 10"^ cm"'^ for isobaric calcula- 
tions started from T = 10* K. Such density, which is relevant 
for the intergalactic medium, is quite low for the interstellar 
medium and, moreover, the cooling time of such gas is larger 
than the current age of the Universe. For higher density iso- 
baric cooling rates calculated from T = 10* K have negli- 
gible dependence on number density. Hence, to study the 
transition of level populations of metal fine-structure lines 
from non-equilibrium to equlibrium we need to take lower 
initial temperature. But this temperature should be high 
enough, that initial difference between isochoric and iso- 
baric cooling rates is small and dependence of the evolution 
on initial temperature is negligible. Using these conditions 
we can constrain the gas temp e rature from se veral x 10"* to 
10^ K (|Gnat fc Sternberg|l2007l : IVasilievll20l'^ ). In this case 
the compression factor (the ratio of initial temperature to 
10^ K) is about < 10^ and we can consider number densities 
n ^ 10^'' cm^^ as initial values. 
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log T, K 

Figure 5. The isobaric (thick Hnes, numeric labels) and isochoric 
(thin lines, lablel 'c') cooling rates for metallicities 10~^ Zq (black 
lines), 1 Zq (grey lines). The numeric labels are the logarithm of 
the initial number density. The initial temperature for the isobaric 
calculations equals 5 X 10® K. The density of isochorically cooled 
gas is 0.1 cm~^. 



Also we should mention that we consider isobaric com- 
pression, which is natural for a gas parcel crossed through 
strong shock front. So the gas temperature ahead a shock 
is much lower than the value T = 3mp«^/16fc, where w is a 
shock front velocity. To mimic this we assume that the ionic 
composition in a gas before a shock corresponds to that at 
T = 2 X 10'* K. Since here we consider pure coUisional ion- 
ization model, we should constrain shock wave velocity by a 
value that does not lead to the radiative precursor formation. 
A stable photoionization precursor wi ll be formed for shock 
veloc ity higher than ^ 175 km s~^ jPopita fc Sutherland! 
1 19961 ) or shock temperature T. = 3mpV^/16fe > 7 x lO'^ K. 
We have found that the cooling rate of a gas behind a shock 
wave with Ts — 5 x 10® K tends rapid ly to the rate in a gas 
coUisionally cooled from T — 10* K (|Vasilievll20l3 ). Thus, 
taking into account the above-mentioned constrains on gas 
temperature we can choose T = 5 x 10® K as initial temper- 
ature. 

In general, the transition from isobaric to isochoric cool- 
ing is controlled by dynamical time of a system: cooling is 
isobaric when the cooling time is much greater then the dy- 
namical time of a system, tc/td 2> 1 and it is isochoric when 
the time ratio is opposite, tc/td <C 1. For a example, the 
evolution of a density fluctuation (e.g., a cloud) is deter- 
mine d by the transition from isobaric to isochoric cooling 
(e.g. , iGnat fc Sternberdl2007D . Such transition is able to de- 
termine the characteristic siz e and mass o f the cold dense 
clumps in two-phase medium jBurkert fc Lin 2000. ). 

Figure [5] presents the isobaric cooling rates for gas with 
metallicity 10"'^ Zq and 1 Zq (thick lines) and the initial 
number densities n = 10^*, 10^^, 10^^, 0.1 and 1 cm^'' (the 
numeric labels in the Figure are the logarithm of the initial 
density). As it is mentioned above the initial temperature 
for the isobaric calculations equals 5 x 10® K. First of all, 
the tails at T ~ 5 x 10® K are co nnected with t he relax- 
ation of gas behind a shock front jVasilievllioi^ ). During 
this relaxation gas is rapidly ionized, the ionic composition 
corresponded initially to T = 2 x 10* K tends rapidly to that 
at r ~ 5x 10'' K. The duration of the period corresponded to 




Figure 6. The electron Xe (upper panel) and molecular hydro- 
gen (lower panel) fractions for metallicities 10~^ Zq (black 
lines), 1 Zq (grey lines) in isobarically cooling gas. The numeric 
labels are the logarithm of the initial number density. The initial 
temperature equals 5 X 10® K. The label 'c' corresponds to the 
isochorically cooled gas with n = 0.1 cm~^. 



almost constant temperature is very short ~ O.Olt/, where 
tf is the total evolution time, which is elapsed from the be- 
ginning to that at which the temperature reaches 10 K. But 
more or less full relaxation reaches during t ~ O.lt/. 

In Figure [5] one can see that the difference between iso- 
choric and isobaric cooling rates reveals at T ^ 10* K. Such 
difference can be explained by two reasons. At first, the re- 
combination of hydrogen and metals becomes more efficient 
in a gas with higher density. At second, the cooling below 
10* K is mainly due to molecule roto-vibrational and metal 
fine-structure transitions, whose level populations go from 
non-equilibrium (L ~ n'^) to the equilibrium {L ~ n) in gas 
with n ~ 1 — 10* cm"''. The first (the hydrogen recombi- 
nation) is significant at T ~ 10'' — 10* K, the second (the 
equilibration) plays a role at T ^ 400 K. 

One can see that for isobarically evolved gas the ther- 
mal instability can develop at T ;S 10* K. The instability 
criterion is w eaker tha n that for isochorically cooling gas 
(|Shchekinov |[l978: Chevalier fc Imamura|[l982l) 



d In T 



(5) 



For Z ^ 10~ Zq this criterion is fulfilled around a small 
bump in the cooling rate at T ^ (4 — 5) x 10'' K, which is 
produced due to the contribution from molecular hydrogen. 
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Figure 7. The contributions to the total isobaric cooling rate 
from major chemical elements for initial number density 1 cm~'^ 
(upper panel) and 10"'^ cm~^ (lower panel) for a gas with Z = 
Zq. The initial temperature equals 5 X 10^ K. 



For gas with Z ^ Zq and initial density n = 10" 
T = 5 X 10^ K the temperature range, where the thermal 
instability can develop, is T ~ (0.2-10) x 10^ K. This range 
becomes narrower with increase of initial density, which is 
explained by steeping of cooling curve (Figure [Sja) due to 
both more efficient hydrogen recombination and equilibra- 
tion of the fine-structure level populations. 

Figure[6]shows the electron (upper panel) and molecular 
hydrogen (lower panel) fractions for metallicities 10"'^ Zq 
(black lines), 1 Zq (grey lines) in gas cooled isobarically. 
The numeric labels are the logarithm of the initial density. 
The increase of density intensifies both hydrogen recombi- 
nation and H2 molecule formation. For solar metallicity gas 
the decrease of ionization fraction compared to that for gas 
cooled isochorically is significant and reaches 0.5-2 order of 
magnitude at T ^ 10'' K (compare the thin dash line to the 
thick ones in the upper panel). So that the higher initial 
density is the larger deviation from the isochoric case can 
be found. The H2 fraction ai Z = Zq increases about 10 
times (thick grey lines in the lower panel), but the H2 con- 
tribution to the cooling rate still remains negligible for solar 
metallicity. We should note that although the H2 fraction is 
higher than that in the isochoric case, but this is still insuf- 
ficient for CO/H2O formation. For 10~^ Zq the dependence 
of ionization kinetics on initial gas density is less pronounced 
(thick black lines in the upper panel). It can be appreciable 
only at T ^ 100 K. The H2 fraction reaches its 'universal 
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Figure 8. The CI and CII fractions in isochorically (solid line) 
and isobarically (dot and dash-dotted lines) cooled gas with solar 
metallicity. In the isochoric calculation the density of a gas is 
0.1 cm-3 (the pressure equals 1.38 X 10-^^{T/W^K) erg cm'^). 
In the isobaric models the initial density values are 10~^ cm~'^ 
(the pressure is 6.9 X 10"^'^ erg cm~'^) and 1 cm"-^ (the pressure 
is 6.9 X 10~^^ erg cm~'^), depicted by dot and dash-dotted lines, 
respectively. The initial temperature for the isobaric calculations 
equals 5 X 10^ K. 



value' and does not depend on density (thick black lines in 
the lower panel). 

Figure [7] shows the contributions to the total isobaric 
cooling rate from major chemical elements for initial num- 
ber densities n = 1 cm^^ (upper panel) and n = 10^^ cm"'' 
(lower panel) for gas with Z = Zq (the contributions by 
N, Ne, Mg and molecules are minor in this temperature 
range for gas with such metallicity). We present the low- 
temperature part of the cooling curve, because the isobaric 
effects reveal in this range. One can note a strong differ- 
ence between cooling rates at T ^ 300 K. The contribution 
by carbon, which is the main coolant for the initial density 
10"^ cm"^, is significantly suppressed for higher initial den- 
sity value. One can find that at T 300 K the energy losses 
in the fine-structure transition of CI (609 iJ,m) dominate in 
the cooling, and also that the critical density (when the level 
populations go from non-equilibrium, L ~ n'^, to equilib- 
L - n) of this trans ition is 160(r/100)"°-^'' cm"^ 



(|Hollenbach fc McKed 1 19891 ). During isobaric contraction 
from 5 X lO'^ K to < 300 K the density increases more than 
> 2.3 X 10^ times. So that at T < 300 K the cooling of gas 
with initial density n — 10~^ cm~^ should be more effi- 
cient (where L ~ than that of gas with initial density 
n — 1 cm""' (where L ^ n). 

As an example we consider the carbon ionization kinet- 
ics in detail. Figure |8] present the carbon ionization states 
for a gas with solar metallicity cooled isochorically (solid 
line) and isobarically (dot and dash-dotted lines). In the 
isobaric models the CI fraction is about 1 at T ;S 300 K. 
Then, there are no other coolants except CI in this tempera- 
ture range (Figure[7]), and the difference between the cooling 
rates (Figure [Tjl is explained by the equilibration of the CI 
level populations. Similar effect can be found for CII within 
T ~ 300 — lO** K, but it is less pronounced than that for 
CI. In general, a decrease of pressure (or initial density with 
fixed initial temperature) leads to that the isobaric cooling 
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rate tend to the isochoric one (see e.g. Figure [5]). For low- 
metallicity gas the role of metals to cooling rate diminishes 
and the transition becomes less distinguishable, because the 
critical density for the H2 roto- vibrational levels is higher, it 
is about 10* cm~^. In the presence of UV radiation this tran- 
sition can be sig nificant in d iffuse clouds in the interstellar 
medium (see e.g. lLisz^l201ll ). 



5 APPLICATION: SUPERNOVA EXPLOSION 

In the previous sections we have considered the cooling rates 
of gas cooled from T 5 x 10^ K. These conditions naturally 
correspond to a gas behind strong shock waves originated in 
collisions of galaxies, virial shocks in clusters of galaxies and, 
more common and usual process in the interstellar gas, su- 
pernova (SN) explosions. The temperature of gas just behind 
SN Shockwave is generally higher than 5 x 10^ K, so that the 
cooling rates calculated in this paper are valid for physical 
conditions behind SN Shockwave. To study the SN evolu- 
tion we use the mu lti-dimensional parallel code ZEUS-MP 
jHaves et al. II2OO6I '). We have added the cooling processes 
into the energy equation using the explicit four-step Runge- 
Kutta method. 

The upper panel of Figure [T] shows that in the high- 
temperature regime (T ^ 10* K) the cooling rate calculated 
in this work is close to that obtained in other studies. A 
significant difference of the rates can be found in the low- 
temperature regime (T ^ 10* K). Certainly, the difference 
between cooling rates more than 2-3 times can lead to ear- 
lier (or later) beginning of the radiative phase of a SN shell 
and, as a result, different scales of fragments appeared due 
to thermal instability. In general, gasdynamic simulations 
with self-consistent chemical kinetics undoubtedly trace any 
variation of cooling rate due to both chemical transfor- 
mation and gas dynamics, but such calculations are very 
time-consuming, restricted in resolution (one should include 
multi-species advection in gas dynamics) and in applications 
(one should take into account valid chemical species and pro- 
cesses for a given physical task). 

Use of tabulated cooling rates gives more flexibility in 



als in the interstellar /intergalactic medium (IFerrara et al. 


20QC 


: Ide Avillez & Mac Low 


2OO2I: iDedikov & Shchekinov 


2004 


: Scalo & Elmegreenl 


2004: Vasiliev et al. 1 120091: 


Ntormousi & BurkertI 201 ll: 


Vasiliev et al. 20121). Metals 



are produced in stellar interiors and ejected into surround- 
ings b y supernova explosions. But available data for cooling 
rates jSutherland fc popitalll993l: ICriat fc Sternberd [ioOTt 



ISchure et al. I l2009l : IWiersma et al. I 120091 ') are very con- 
strained in metallicity and temperature, o r even combined 
from different sources jSchure et al. Il2009^ . Of course, one 
can utilize pre-compute d tables assumed the equilibrium 
(e.g. use CLOUDY, see IWiersma et al. I |2009| '1 or choose 
several species for nonequilibrium calculati ons, whereas 
other s are considered in equilibrium (e.g. ISmith et al. I 
|2008| ). But in the former we have equlibrium rates, which 
significantly differ from the nonequilibrium ones, whereas 
the latter needs a significant computational time. 

Besides our cooling rate s Figure [T] presents two other 
cooling curves at T < 10* K : ISpaans fc NormanI (|l997l ') and 
iDalgarno fc McCravl (|l972h . The former differs slightly from 



the rate calculated here at T ^ 10^ K, so we don't con- 
sider it. The latter widely used in the ISM studies (e .g. , 
Ide Avillez fc Mac Lowll2002l : iNtormousi fc Burkertll201ll ) is 
valid only at T < 10* K and demonstrates significant dif- 
ferences from ours for Xe ^ 0.1. So that we use the cooling 
rates calculated in this paper in high temperature range, 
T > 10* K, whereas for T < 10* K we take three options 
for cooling rate: a) ou r self-consistent isochoric ca lculations 
(label 'NEQ'), b) the lOalgarno fc McCravl (| 19721 ) rate for 
the fixed electron fraction Xe = 0.1 (label 'NEQ+DMC72, 
Xe = 0.1') and c) the same rate, but for Xe = 10~* (la- 
bel 'NEQ-(-DMC72, Xe = 10"*'). Our choice of x^ is ex- 
plained that the rate for higher value Xe is close to our rate 
at T ^ 10^-10* K , wher eas the lower Xe is close to the value 
that lSmith et al~l (j2008l ) found for the interstellar gas with 
solar abundance (the same value can be obtai ned for gas ion- 
ized by cosmic rays in the local ISM, see e.g. IWolfire et al. I 
Il995h . 

Here we investigate the dependence of the SN shell evo- 
lution in a homogeneous medium with solar metallicity on 
different cooling rates in low temperature range, T < 10* K. 
We have carried out the 2D gasdynamic simulations (cylin- 
drical geometry) of a SN explosion. The numerical resolu- 
tion is 1000 X 1000, that corresponds to the physical size 
of a cell 0.4 pc. The SN energy is initially injected in the 
central region with radial size 3.2 pc, the total energy is 
lO"^^ erg. To exclude the evolution of the gas in the medium 
before Shockwave we set the temperature of this gas to 50 K 
and do not allow to cool it down. The differences of the SN 
shell evolution due to variation of the cooling rate can be 
found just after time t ~ lO'^^ s, when initially hot shell 
cooled to ^ 10* K. We continued the calculations as long 
as t = 10^* s, at this time various instabilities destroy SN 
shell, but here we are interested in the beginning of the shell 
evolution cooled below 10* K. At i ;S 3 x 10^^ s the shell is 
still close to spherical shape, so that we can present radial 
distributions of physical quantities only. 

Figure |9] present the radial distributions of density (up- 
per) and temperature (lower) around a SN exploded in a 
homogeneous medium with density n = 3 x 10~^ cm~^ and 
solar metallicity for different cooling rates at T < 10* K. 
At t = 1.4 X 10" s (the left row of panels in Figure [SJ 
one can find that the radiative phase begins in the SN 
shell co oled with the rate ca l culate d in this work ('NEQ') 
and the iDalgarno fc McCravl (| 19721 ) rate for the fixed elec- 
tron fraction Xe = 0.1 ('NEQ-I-DMC72, Xe = 0.1'). In the 
first ('NEQ') model the thin shell is formed, the density in 
the shell increases in several times and the temperature is 
low as ~ 100 K. Whereas in the second ('NEQ-I-DMC72, 
Xe = 0.1') model the effects of radiative cooling begin to 
become apparent only: the radiative phase have just be- 
gun. The main difference between two models is the absence 
of temperature plateau at T ~ 10* K in case of the self- 
consistent cooling rate ('NEQ'). The plateau (see radial dis- 
tances r ~ 200—250 pc) originates from lower cooling rate in 
the 'NEQ-I-DMC72, Xe = 0.1' model in comparison with the 
rate used in the 'NEQ' model. At t = 1.6 x 10^^ s the SN shell 
becomes radiative for the 'NEQ-I-DMC72, Xe = 0.1' model. 
A gas in the shell begins to cool below 10* K, and the plateau 
in the radial temperature distribution disappears gradually. 
In the third ('NEQH-DMC72, Xe = 10"*') model considered 
here a gas in the shell starts to cool down below 10* K only 
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Figure 9. The radial distributions of density (upper panels), temperature (lower panels) around a SN exploded in a homogeneous 
medium with n = 3 X 10~^ cm~'^ and solar metallicit y for cooling rates in low t emperature range, T < 10* K: a) the self-consistent 
calculations presented in this paper (solid lines), b) the lDalgarno fc McCravl |l97j rate for the fixed electron fraction = 0.1 (dashed 
lines) and c) the same rate, but for Xe = 10~* (dotted lines). Left column panels shows the distributions at t = 1.4 X 10^'' s after 
explosion, middle column presents the same at t = 1.6 X 10^^ s and right column of panels shows the same at t = 2.4 X 10^'^ s. In high 
temperature range, T > 10* K, we use the cooling rates calculated in this work. 



at t = 2.4 X 10^^ s (right column of panels in Figure[9|. This 
delay can be easily explained by lack of cooling at T ~ 10* K 
in the 'NEQ-I-DMC72, Xe = 10"*' model in comparison with 
the other models considered in this section (Figure [!}. Note 
that this delay can significantly affect on the SN hot bub- 
ble dynamics: the hot bubble has different sizes depending 
on the cooling curve assumed in calculation. This can be 
apparent in mixing metals, because the metals ejected by 
SN are initially confined in the hot bubble. Thus, we can 
conclude that tabulated cooling rates should be used care- 
fully, because some gaseous structures may form due to using 
cooling rates cooling rates calculated for different physical 
conditions. As far as possible one should use cooling rates 
calculated self-consistently for the whole temprature range 
of one's interest. 



6 CONCLUSIONS 

In this paper we have presented self-consistent calcu- 
lations of the non-equilibrium cooling rates of a dust-free 
collisionally controlled gas with metallicities in the range 
10~* Zq ^ Z ^ 2 Zq cooling from high temperatures 
T > 5 X 10^ K down to 10 K0. We have found that 

• molecular hydrogen dominates cooling at 
10^<r<10* K and Z < 10~^ Zq, its contribution 



^ The isochoric cooling rates presented in Figure[2j3 are available 
in electronic form at |http:/ /ism.rsu.ru/cool/coolfun.html[ 



around T ~ (4 — 5) x 10^ K stimulates thermal instability 
at Z < 10"^ Zq; 

• the ionization fraction remains rather high at 
T < 10* K: it reaches up to ~ 0.01 at T ~ 10^ K and 
Z 0.1 Zq, and becomes higher for higher Z; 

• despite high ionization fraction the abundance of molec- 
ular hydrogen decreases with metallicity; 

• at T ^ 10* K isobaric cooling rates are lower than iso- 
choric ones; 

• in isobaric processes the ionization fraction decreases 
considerably and molecular hydrogen grows more efficient 
compare to isochoric ones. 

We have compared our self-consistent cooling functions with 
those published in previous works and found considerable 
differences. In particular, we found that the two approaches 
for the cooling processes result in a fairly distinct dynamics 
of supernova remnants. 
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APPENDIX A: MOLECULAR HYDROGEN 
AND DEUTERIUM CHEMISTRY 

In this section we present the minimal model of chemical 
kinetics for molecular hydrogen and deuterium in the low 
temperature range, T < 2x 10*K. We include colUsional gas- 
phase reactions only. The reactions are Usted in Table EH 
This model is based on the well-tested minimum models for 
molecula r hydrogen and d e uterium chemistry for both pri- 
mordial (lAbel et al. Ill997l: ICalU & Pallalll998l) and metal- 
enriched ( Omukai et al. II2005I : Glover & Jappsenll2007^ gas. 
We don't list here the helium chemistry, because it consists 
of coUisional ionization and recombination only, these pro- 
cesses are inevitable part of our model and the references 
to reactions can be found in (|Vasilievll20"lll ). Other helium 
species, like HeH^, cannot change thermal evolution of a gas 
significantly, so that we do not include its in the model. 
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Table Al. Chemical reaction rates 



reaction Reference 



H2 chemisrty 


H + e- H- 4 


■ hv 


GaUi & Palla (1998) 


H- + hi^ ^ H + 


c~ 




H- + H^H2 4 


- c" 




H- + H+ -> H+ 


4e- 




+ H+ ^ 2H 






H+ +H ^. H+ - 


h \iU 




H+ + H ^ Ha 4 


-H+ 




H+ 4e- ^ 2H 






H2 4 H+ ^ H+ 


4H 




H2 4 e- ^ 2H - 


he- 




H2 4 H 3H 




Mac Low & ShuU (1986) 


H+ 4c- ^ H4 


hv 


Verner et al. (1996) 


H4e- -> H+ 4 


2c- 


Abel et al. (1997) 


H- 4e- ^ H4 


2e- 




H- 4 H ^ 2H4 


- e~ 


ShaDiro & Kane (1987) 


D chemistry 


D+ 4c- ^ D 4hi/ 


GaUi & Palla (1998) 


D 4 hi/ ^ D+ 4 


e 




D 4H+ D+ 


l-H 




D+ 4H D 4 


H+ 




D4H2 ^ H4HD 




D+ 4 H2 HD 


4H+ 




HD 4 H ^ H2 4 


- D 




H+ 4 HD H2 


4D+ 




D 4e- D" 4 


- hv 




D+ 4D- -> 2D 






H+ 4D- D- 


l-H 




H- 4 D ^ H4 


D- 




D- 4H^D4 


H- 




D- 4H^HD- 


fc- 
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